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A new formulation is presented to treat thermal and quantum fluctuations around the 
kaon condensate on the basis of chiral symmetry; separating the zero-mode from the be- 
ginning and using the path-integral method, we can formulate the inclusion of fluctuations 
in a transparent way. Nucleons as well as kaons are treated in a self-consistent way to the 
one- loop order. The effects of the Goldstone mode, stemming from the breakdown of V — 
spin symmetry in the condensed phase, are figured out. 
, A procedure is discussed to renormalize the divergent integrals properly up to the one- 

■ loop order. Consequently the thermodynamic potential is derived. It is pointed out that 

y—t ' the zero-point fluctuation by nucleons gives a sizable effect, different from the kaonic one. 
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I. INTRODUCTION 

> : 

It has been extensively discussed that kaons, the lightest hadrons with strangeness, show characteristic 
^^"O features through the kaon-nucleon interactions once immersed in nuclear medium. The scalar interaction 
(called the sigma term) and the vector interaction (called the Tomozawa- Weinberg term) are the leading 
ones, and they bring about a large attraction for the excitation energy of K~ mesons at low momentum, 
whereas a repulsion for K + mesons. There are many proposals to see these consequences experimentally; 
, K~ enhancement in the subthreshold kaon production or the anomaly around the meson peak in the 
dilepton production from the relativistic heavy-ion collisions are some of them [|-|. 

In 1986 Kaplan and Nelson suggested a possibility of kaon condensation in dense matter by using a chiral 
Lagrangian Q|. Since then this subject has been studied by many authors j^]. The s-wave attractions 
between kaons and nucleons remarkably reduce the excitation energy of K~ mesons in neutron matter and 
thereby kaons should condense when the degeneracy energy of electrons exceeds it. It is well-known that 
the kaon condensation results in the large softening of equation of state (EOS) and the inducement of new 
processes due to the presence of the condensate. Since it may occur at relatively low densities, (3 — 4)/?o with 
Po being the nuclear density ~ 0.16fm -3 in neutron star matter, it is believed to have many implications 
on neutron star physics; the fast cooling mechanism, rapid rotation or the low-mass black hole scenario are 
some of them [0-p). 

We have studied this subject from the point of view of chiral rotation on the chiral manifold G/H with 
G = SU(3)l x SU(3)n,H = SU(3)v ||- We have defined the kaon condensed state by acting a unitary 
operator U on the normal matter |0) (the meson vacuum); \K) = Uk((0))\0) with Uk((Q)) — exp(i(0)F| ), 
where F| is the 4-th component of the axial-vector charge of the SU(3)l X SU(3)r algebra and (8) is the 
chiral angle. [] Since the kaon condensed state is a chirally rotated one, the change of the kaon-nucleon 
dynamics from that in normal matter arises from the symmetry breaking terms in the Hamiltonian. Thus 
we have succeeded in extracting the essential features of kaon condensation in this formalism. 

So far almost all the studies have been done within the mean-field theory, the classical approximation. 
Recently there begin to appear studies about kaon condensation in protoneutron stars, where thermal effects 
or roles of the neutrino trapping become important and we must treat it beyond the mean- field level @. 
Low-mass black hole, if exists, may be produced in this era |8j]. A numerical simulation based on the general 
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relativity has been already done about the delayed collapse of protoneutron stars to the low-mass black hole 
|PH . However, they used a cold EOS for the kaon condensed phase, while temperature (T) rises to several 
tens MeV there. Also, in this era, various time-scales become important; collapsing, neutrino-trapping or 
initial cooling time-scale may be relevant. Hence we cannot assume chemical equilibrium a priori for such 
hot matter since it may be achieved through weak interactions. A problem how condensate appears and 
grows to establish chemical equilibrium is an interesting subject to pursure fl(][| . Prakash et al. treated the 
kaon-condensed phase at finite temperature and discussed the properties of protoneutron stars within the 
meson exchange model, since there is no consistent theory based on chiral symmetry |T^| . 

In a recent paper, Thorsson and Ellis have tried to include quantum or thermal fluctuations within the 
heavy-baryon chiral perturbation theory (HBCPT) fl5|, frozing dynamical degrees of freedom of nucleons to 



the thermodynamic potential 
of kaons at T = to the one 



13fl . However, they have discussed only the effect of the zero-point fluctuation 
oop order. The most interesting result in their paper may be the dispersion 
relation of kaons in the condensed phase, which is a key ingredient for the discussion of thermal properties 
of the condensed phase. Moreover, it plays an important role in the discussion of the relaxation process 
from the nonequilibrium state to the kaon condensed state in equilibrium fiof . They have used the Kaplan- 
Nelson Lagrangian, where the basic variable is the SU(3) matrix specified by the eight Goldstone fields 
4> a , U(4> a ) — exp(2iT a (f> a ) with T a being the generators of SU(3) Lie algebra. They picked up only the 
degrees of freedom for the charged kaon fields K , namely 04 and 05, and separated the classical field (the 
condensate) and the fluctuation field around it in the standard manner, K = (K ) + K . The resultant 
dispersion relation for kaons has a very complicated form due to the interaction between the condensate 
and fluctuations. Accordingly other thermodynamic quantities also become too complicated to be tractable. 
Thus there is no consistent finite temperature calculation for the chiral case. 

In ref. fll4| we have presented a new formalism to treat the kaon condensation at finite temperature within 
HBCPT. We, in this paper, extend the formalism to take into account the dynamical degrees of nucleons as 
well as kaons |l(|. We introduced fluctuations around the condensate by acting the chiral transformation 
on the chiral manifold successively; after generating fluctuations around the meson vacuum by a chiral 
transformation (77) we transform this state further by another chiral rotation (£) , which should be the same 
as the classical one (see Fig. I). 




FIG. 1. Chiral rotations on the SU(3) manifold. Circles C, CK indicate the fluctuation "areas" around the vacuum 
(O) and kaon condensed state (K), respectively. 



Since the constant condensate can be considered as a kind of zero-mode, this procedure can be regarded 
as a separation of the zero-mode from the original matrix U{<j) a ). Similar procedure has been applied in the 
different context to see the finite- volume effects of the chiral symmetry restoration . Then we can obtain 
the dispersion relation for the kaon excitation around the condensate. We have seen that the dispersion 
relation shows the similar feature to the one by Thorsson and Ellis, whereas its form is very different. 
We have also proposed an approximation which makes the expressions of the thermodynamic quantities 
considerably simple. We have checked numerically that this approximation may work very well. 

In this paper we treat the dynamical degrees of freedom of kaons and nucleons consistently; since nucleons 
interact with kaons, their single-particle energies are modified not only by the condensate but also by the 
fluctuations. For this purpose we use the Hartree approximation for the four-point vertices among kaons 
and nucleons; the nucleon Green function contains kaon loops and the kaon Green function does nucleon 
loops. Therefore we sum up the infinite series of bubble diagrams composed of one-loop diagrams by kaons 
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and nucleons in a self-consistent way. We shall see that our formalism makes the analysis of the structure 
of kaon-nucleon (KN) dynamics transparent and clarifies physics included there. 

Generally these loop diagrams include divergences, which should exist even at T — 0. We discuss how 
finite contributions can be extracted in our formalism. Analyzing the structure of the self-energy terms, we 
show these divergences are properly renormalized to the one-loop order by introducing relevant counterterms 
corresponding to the nucleon and kaon masses, and the KN and KK interaction vertices. Then no more 
counterterms are needed to renormalize the thermodynamic potential. Effects of the zero-point (quantum) 
fluctuation by kaons has been shown to be small (T^JIJ]. We shall see that the zero-point fluctuation by 
nucleons has a sizable effect on the ground-state property of the condensed phase. 

In H we present our formalism to treat fluctuations around the condensate by way of the path integral 



approach. In [II we develop our approach and perform the one-loop calculation. There we get the dispersion 



relation for kaons which is a key object to discuss the properties of the kaon condensed phase at finite 



temperature. We discuss in IV how the thermodynamic potential is renormalized by analyzingthe self-energy 
terms. We also give some numerical results about the zero-point fluctuation by nucleons. |y| is devoted to 
some discussions and concluding remarks. In Appendix A we discuss the separation of zero-modes, one of 
which is the would-be condensate. Some results within the mean-field approximation are summarized in 
Appendix B. In Appendix C we give some properties of the kaon Green function and expressions of the loop 
contributions. 



II. PATH INTEGRAL FORMULATION 
A. Partition function 

We start with the Kaplan- Nelson Lagrangian as a chiral Lagrangian 0] , 

C c hiral = C + C$B, (2-1) 

where Co is the symmetric part, 

C Q = ^tr(^ I C/ t ^C/)+ztr{B(/?-77i B )S} + mr(B7 M 7 5 {^,S})+^tr(B 7Al75 K,B]). (2.2) 

and Csb the symmetry breaking part; 

C SB = vtvm q (U + -2) 

+ aitr£?(£m g £ + h.c.)B + a 2 trBB(£m q £ + h.c.) + a 3 {trBB}tr(m q U + h.c.) (2.3) 

with the quark mass-matrix, rh q ~ diag(0, 0, m s ). The coefficient v implies the vacuum expectation value 
of the quark bilinear, — v = (0\qq\0), and the coefficients measure the strength of the explicitly symmetry 

breaking: the kaon mass is given as m 2 K ~ vm s / f 2 and the KN sigma terms as E/fp = —m s (ai + + 203) 
and Sjfn — ~^s(fl2 + 203). U is the S£/(3) matrix parametrized by the eight Goldstone fields (f> a , 

U = cxp[2zT a a //] e G/H (2.4) 

and Vfj,, are the vector and axial- vector fields constructed by U, 

V^l/2(£t^ + ^t) 

Ap = i/2(Z%t (2.5) 
with U — £ 2 . The baryon octet B can be represented by the ST/ (3) matrix, 
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Then the covariant derivative D„ is defined by the vector field V M ; 
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D^B = d^B + [V^B}. 



(2.6) 



The transformation properties of the fields U and B under G are found in, e.g., ref. p4| . 

To treat thermal and quantum fluctuations consistently we use here the path-integral formulation. It is 
well-known that chemical potentials /i a , which should be taken into account to ensure various conservation 
laws in the ground state, can be introduced by "gauging" the Lagrangian with artificial "gauge fields" = 
T a A£ with A^ = (/i a ,0). Here we must take into account conservation of two quantities: electromagnetic 
charge and baryon number. Accordingly we consider two chemical potentials, the charge chemical potential 
and baryon chemical potential. For the charge chemical potential fj,K, which we shall see to be identified as 
the kaon chemical potential, the "covariant derivative" V^JJ for the Goldstone field then reads 

V li=0 U = dU/dt + i f x K [T em ,U], V^U = d^U for others, (2.7) 

with the charge operator, T em = T3 + 1/V378 = diag(2/3, — 1/3, — 1/3). For baryon field B, we must 
introduce another Abelian gauge field to ensure the conservation of baryon number, B^ — (/i„,0) with the 
baryon chemical potential fi n , which we shall see to be the neutron chemical potential. Then the "covariant 
derivative" W^ = qB for the baryon field reads 

V^B = dB/dt- in n B + i m [T em ,B] + [V ,B], W^B = D^B for others. (2.8) 

Therefore the partition function in the imaginary-time formalism can be represented as follows (r = it, (3 — 

l/T); 

Zchirai — N f [dU] [dB] [dB] exp[S' chiral ] (2.9) 



S' chiral = dr J d 3 xC chlra i (V„U, U, V„B, B) = ^ ' dr J d 3 x[C chlra i + SC], 



(2.10) 



where the "time-derivative" should be read as the derivative with respect to imaginary time. Here we can see 
that there appears an additional term SC besides the original chiral Lag rangian £ c hirai , which is in general 
non-invariant under the chiral transformation. For the Lagrangian (2.1), we find 

SC = -^tr{[T em , U}^-+ f^[T em , C/t] } _ ^£ tr{B t[ (e t[ TemjC ] + £[T em ,£t]), B }} 
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L^K t r{[T em ,U}[T em ,U^]} + fi n tr{B^B} - m ti{B^[T emi B]}. (2.11) 



Here the first term corresponds to the mesonic charge, and the last two terms indicate nothing else but the 
baryon number and the electromagnetic charge of baryons. 



B. Transformation of the coordinates 



We show our formulation to include thermal and quantum fluctuations in the chiral Lagrangian. Before 
that let us recall the tree case. We have described the kaon condensed state \K) , which is specified by the 
chiral angle (9), the order parameter of kaon condensation, as a chiral- rotated state ||; 

\K) = U K ((0))\O) (2.12) 

with the operator Uk{(9)) = exp(iF^ (9)), where Ff is the 4-th element of the axial-vector charge in 
SU(3)l x SU(3)r algebra and |0) represents the normal nuclear matter. We can also describe this 
point on the SU(3)l x SU(3)r/SU(3)v — SU(3) manifold which is coordinated by the Goldstone fields, 
U(4> a ) — exp(2iT a a //) £ SU (3) and the meson vacuum |0) corresponds to the fiducial point on the manifold 



by the spontaneous symmetry breaking. Then kaon condensed state Eq. (2.12) corresponds to the SU(3) 
matrix Uk = exp[2i{T4(</>4) + 25(^5)}//], and can be represented as a chiral rotated one from the vacuum 
on this manifold, 

U K {{9)) = (UyC = C 2 (2.13) 
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where Uy = U(<f> a = 0) = 1 corresponds to the meson vacuum and £ is the operator, ( — exp(i(M)/V2/) 
with the constant matrix 



M = 



0^+ 

K~ 



K* = (d>4±io 



and 9 2 = 2K+K-/f 



Thus we have seen that the classical kaon field (K ± ) induces the chiral rotation on the manifold 



l=exp(2iT a ' 0/f) 
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U f = n . T| =exp(2iTa?a / 
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:FK 



FIG. 2. Diagram showing the relation among the transformations on the chiral manifold (c.f. Fig. 1). The symbol 
TE indicates the one by Thorsson and Ellis [13]. 

When we introduce the quantum or thermal fluctuation, there may be several ways (Figs. 1,2). Thorsson 
and Ellis introduced the kaon fluctuation field K as a deviation from the classical kaon field in the standard 



manner, 



(2.14) 



Here we introduce the fluctuation fields by developping the idea of chiral rotation mentioned above. First, 
we notice that any fluctuation of the Goldstone fields can be introduced by a chiral rotation Ut from the 
vacuum Uy = 1, 



■q : Uf = -qUyrj = rf 
with rj — exp(iT a <p a / /). The subsequent chiral rotation by £ transforms Uf into the form, 

(: [/«0),0 Q )=CWa)C- 



(2.15) 



(2.16) 



We can easily see that in t he lim it <j) a —> 0, which corresponds to the previous classical approximation, the 
matrix U is reduced to Uk ( |2.1S ). Thus we have introduced the fluctuation fields by the two-step procedure. 
^ It is to be noted that the constant component in <fi a is redundant since it can be always absorbed into £ 
by redefinition. It may be also worth noting that we can easily take into account fluctuations of any meson 
field in the kaon condensed state by the form of Eq. ( 2.1 6| ) . Q Accordingly, £ operator, which is defined by 
U = £ 2 , can be obtained by solving the subsidiary equation, 



(2.17) 



where the matrix u is defined by the second equality in Eq.(2.17). It depends on <j> n , (9) in a complicated 
nonlinear way, and thereby it is very difficult to find u for the general form of Uf J24J. A systematic way 
to find the form of u may be the perturbative method: expanding the matrices 77, u with respect to the 



2 Hereafter, we omit the tilde on the field for the sake of simplicity. 

3 The similar decomposition was also used in the context of the loop expansion within chiral perturbation theory 
[H or the bound-state approach to the Skyrme model 

Alternatively, we can regard this procedure as a separation of the "zero-mode" from the full SU(3) matrix U(4> a ) 
(see Appendix A). Using Eq. ( |A17j ) we can see that the value of the thermal average (9) should be determined by 
imposing the extremum condition for the resultant thermodynamic potential at the end. 
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fluctuation field <f> a , we can solve Eq. (2.17) order by order. Hence we can find the relevant form of u up to 
any order of <p a , depending on the problem or the approximation we are interested in. On the other hand, 
if we restrict ourselves to the fluctuations of kaon sector s.t. 

U f = exp(iv/2M//), (2.18) 

then we can find u in the closed form, 

u = diag(K*/|K|,l,«/|*s|), (2.19) 

with 

K = cos«0>/2) cos(0/2) - j^^j sin«0)/2) sin(0/2). (2.20) 
Here the following relation is to be noted, 

u^T ern u — T em . (2.21) 

Defining a new baryon field B' by the use of the matrix it, 

B' = u^Bu, (2.22) 

we can see that the chiral-invariant pieces of the Lagrangian are not changed and only the symmetry breaking 
pieces play essential roles, 

C chira i(U K ,B)=Co{U K ,B)+CsB(U K ,B) — > C (U f , B') + C SB {(U f (, uB'vt) 

5C(U K ,B) — ► SCiCUfCuB'vt). (2.23) 

Thus all the dynamics of kaons and baryons in the condensed phase are completely prescribed by the non- 
invariant terms under chiral transformation. It is to be noted that this feature is quite similar to the one 
within the classical approximation at T — Q. 
Using Eqs. ( pl9| ), p!2^ ), we find 



6£((U f (,uB'ui) = -(cos«0) - l)fi K tr{B^[V 3 ,B'}} + f 2 f i 2 K /2 S m 2 (9) 

- hk cos(9)(K + K- -K+K-) +fi K cos{6)K + K-/ftr{B'i[V 3 ,B']} 

+ cos 2 (9) fi 2 K K + K- - fi 2 K / 4 sin 2 (9) (cxp(i"f)K- + exp(~i-/)K + ) 2 

+ M„tr{S'tS'} - » K tr{B'i[T em ,B'}}..., (2.24) 

with the V-spin operator, V3 = 1/2(T3 + \/3Tg), and K ± = dK^/dr, where the phase 7 is defined by 
exp(±«7) = (K ± ) /\(K)\, which can be absorbed into the fluctuation field by redefining it, exp^i^K^ 1 — > 
K . In other words the phase of the condensate is arbitrary, reflecting V3 symmetry in the chiral Lagrangian, 
and we can choose 7 = without loss of generality. As a consequence V3 symmetry is spontaneously broken 
in the condensed phase, and thereby we shall see the Goldstone mode associated with this symmetry (see 

0). Q 

If Csb originates from the quark mass term , C mass = —qrh q q, we can further explore the transformation 
property of Csb in a model-independent way: first we can see 

Cmass -> C m ass ~ t SU1 9{F% , C m ass] + (cOsO- l)Sif (2.25) 



with the sigma term Sj<- = [F% , C mass ]], under the transformation by r] in Eq . (2.15). Similarly, C ma ss 
transforms like (2.25) with (9) instead of 9 under the chiral rotation £ in Eq. ( 2.1 6| ) . Hence C m ass eventually 



transforms under the successive chiral transformations (2.15) and ( 2.1 6| ) like 



Cmass C mass - i (sin((0) + 9) - sin0) [Fl,C m ass] + ((cos((fl) +9)- cosff) E*, (2.26) 



3 We met the similar situation in the context of pion condensation 121 
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where we have used the property [F% , [F% , [F%, C ma ss]]] — [F$ , jC mass ] , The second term is a parity-odd 
operator and irrelevant due to the fact, (975(7) = 0. The term proportional to sin(#) included in the la st ter m 
is also irrelevant in the following calculation bacause of the thermal-equilibrium condition (c.f. Eq. (A17)). 
Therefore there is essentially left the term, 



Cmass -> Cmass + (cOs(6>) - 1) COS 9T, K . 

Accordingly, in terms of the effective Lagrangian £5 b , it should transform like 

Cs B {QUf(;,uB'v}) = jC S B(Uf,B') + (cos<0) - 1) cos0£ K (l, £') 
besides the irrelevant terms, where £#-(1, B') has a standard form, 

£k(1, B') = f 2 m 2 K - Y>k p P'p' - T, Kn n'n', 



(2.27) 



(2.28) 



(2.29) 



for nuclear matter. [] In Eq. ( 2.29| ) we have used the following notations: f 2 m 2 K — v[m u + m s ], £k p („) = 
(p(n)[EiK \p( n ))- This consequence is rather general and, of course, consistent with that in the Kaplan-Nelson 
Lagrangian (2.3). 

For the integration measure (N = 3), it is given as 



N 2 



[dU] = M{cj>) 11 [d<f>i 



(2.30) 



where M{4>) is the Lee- Yang term |||, M(4>) = exp ^ (4) (0) ffidr Jd 3 x\ng((t>) with the Cartan-Killing 

metric on the SU(N) manifold from the invariant line element, ds 2 = gik(<t>)d4>id<j)k = l/2tr[c?L r <i[/^l. It is 
to be not ed th at the Lee- Yang term is important to make the integration measure chiral invariant. f\ The 
measure ( 2.30 ) is invariant under the transformation ( |2.16| ) by £, 



[du] = [d((u f ()) = [du f ], 



(2.31) 



because G is the isometry, tx[C,dU /CC^ dUttf] — tr[dU fdlll], so that the geometry around the condensed point 
is the same as in the vicinity of the vacuum. Since we are concerned with the fluctuations, the measure can 

be approximated by the one with flat curvature ]l9| , [dUf\ ~ Y\a=i 1 [d ( t>i\- I n this case we are not worried 
by the Lee- Yang term because of g = 1 + 0(4> 2 /f I ). Thus we find 



N- 



Wi] [dB'\ [d&] eMS e J h L(C, U f , B') 



(2.32) 



i=l 



III. THERMODYNAMIC POTENTIAL 



A. Partition function for a chiral Lagrangian 



Using the formulation given in |H| we evaluate the partition function in the kaon condensed state. ^ In 
the following calculation we only consider the one-loop diagrams, and thereby we retain only the quadratic 



3 We, hereafter, omit the prime on the nucleon fields for simplicity. 



In ref. p3| , their treatment of the Lee- Yang term is obscure: 7 1 / 2 in Eq (13) in ref. [[[3), which corresponds to g in 



our notation, should read exp 



i<5W(0)/;dr/d 3 xln 7 



We briefly discuss the case of the classical approximation in Appendix B, where we can see that our previous results 
given in ref. M are recovered in a model-independent way. 



7 



terms with respect to the kaon field. Moreover, we are concerned here with only charged-kaon fluctuations 
and nuclear matter without hyperons. [] Then we can write 



(3.1) 



with the spinor ip 



where 



^ ] by the use of Kaplan-Nelson Lagrangian ( fTI 

S chiral = Sc + S K + S N + S int , 



S C =[3V 



-H 2 K f 2 sm 2 (9)-f 2 m 2 K (l-co S (0)), 



(3.2) 



(3.3) 



dr / d x 



d u K + d u R- - cos(6)(m 2 K - fi 2 K cos(9))K+K- - fj, K cos(6)K + d T K~ 



-^f sm 2 (e)(K+ + K-f + 0(\K\ 4 ) 



with the shorthand notation d u K + d u K" = -d T K + d T K~ - VK + VK~ 



S N = I <It j d\r 



d 

+ /Xif (1 - cos(#)) 7 ^-^ + (1 - cos{6))V KN )i> 



(3.4) 



(3.5) 



and Si n t is the interaction term, which mainly consists of the two terms, the sigma-term and the Tomozawa- 
Weinberg (TW) term, 



Sint = j dr I d".r 



2 cos(9}^T, KN ipK-K^ 



1 
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(3.6) 



Hn are 2x2 diagonal matrices in the isospin space: Y>kn 



with the KN sigma term Yikn- It is to be noted that the sigma-term Y^kn and nucleon chemical potential 

Ek„ ) aild ^ = ( O" Hn ) 

fJ-p = Mn - ^Ki respectively. 

Since the effective action S^ ral is already bilinear with respect to the nucleon field, we can easily integrate 
out the nucleon degree of freedom in Eq. (3.1), 



ZcMral = N J[dfa}[dfc] exp[S' K ]DetG = N J [d^ A ][dfc] exp[S^ + TrlnG]. 



(3.7) 



with S' K = S c + Sk and the Dirac operator G. The capital letter for determinant or trace is used for the 
Dirac operator G of infinite dimension over space-imaginary-time. The Dirac operator G is composed of two 
parts, G — Go + SG, with the "free" part 

Go = + ^7o7 ■ V - Mj + m+ A*x(l - cos(0))^-^ + (1 - cos(6))j ^kn, (3.8) 

OT 4 



9 An extension to include hyperons may be straightforward, but note that there is still large ambiguity about their 
coupling strength with hadrons, e.g. see ref. Jl^ ]. 
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and the interaction part 

5G =^fco S (9) 7Q K-K+ - "^^p \k + d T K~ - 2 t x K cos{6)K+ if" J + 0{\K\ 4 ). (3.9) 

A naive treatment may proceed with treating Go as the zeroth-order operator and expanding TrlnG in 
terms of <5G, 

TrlnG = TrlnG + Tt[Gq 1 6G] + 0(\K\ i ). (3.10) 

The first term gives a one-loop contribution of nucleons to the partition function, and the second term to the 
kaon self-energy term. The latter, together with Sk, gives a one-loop contributions of kaons to the partition 
function. However, this procedure is not adequate for treating the KN interactions in a consistent way; a 
simple analysis of the loop counting shows that the kaon one-loop contribution includes an infinite series 
of nucleon loops, while the nucleon one-loop contribution does not include any kaon loop. In the following 
we use the Hartree approximation to take into account all the bubble diagrams by both nucleons and kaons 
non-perturbatively, while the exchange diagrams are discarded. Generally speaking, the exchange diagrams 
are suppressed in the high-density or high-temperature limit, compared with direct diagrams (see Fig. 3 for 
the three loop case) f^Q. 



(a) (b) (c) 

FIG. 3. Three-loop diagrams, (a) and (b) are the direct (bubble) diagrams, which are taken into account by 
the Hartree approximation, (c) is the exchange diagram, which should be suppressed by the power of density or 
temperature in comparison with (a) or (b). 

We define a new operator Gh, Gh = Gq + (5G) by taking the thermal avarage of 5G for kaons, 

G H = + ijol -V + fi N + ^(n^ + 2fi K (l - cos<0») - M* 7 ° (3.11) 

/ M* \ 

with the effective mass of the nucleon, M* — I p , , rsB with 



M* 

M* = M - T, Ki /f cos(0)nk - T, Ki (l - cos(0)). (3.12) 

The quantities n x K ^ K mean the thermal kaon-loop contributions for the self-energy of nucleons and they 
can be represented in terms of the thermal Green function of kaons (see Appendix C) ^] ; n x K = (K + K~) 

and n 2 K = —n° K + 2fi K cos(9)n 1 K with n° K = (K + d T K~). 
Using Gh we can rewrite Eq. ( 3.1C ) as 

TrlnG = TrlnG H + Tr[G ff ^G] - Tt[G h 1 {8G)] + 0{\K\ A ). (3.13) 



10 The quantity v} K includes the divergence and should be properly renormalized. We shall meet other quantities 
that also include some divergences in this section. We leave the renormalization problem in 
them schematically in this section. 



IV 



and hereafter treat 



9 



Substituting Eq. (|3.13|) into Eq. (pO), we find 



Z ch i ra i =N J [d0 4 ][dfo] exp[^]DctG H , (3.14) 

with 5 A - = S c + S e J f - Tr [G" 1 (6G) } + ( | K | 4 ) . Thus we construct an effective action for the kaon fluctuation 
field, 

S e J f = dr J d 3 x [-d^K+d^K- + cos(6»)(^ A - cos(6>) - m 2 K + a + 2b^ K )K+K- 

- (hk cos(0) + 6)^+ 3 T X- - ^ sin 2 (0)(X+ + if") 2 + 0(\K\ 4 ) ] , (3.15) 

where a = f- 2r Tt[G^^T, KN \ and b = /- 2 Tr[G^(3 + r 3 )/8] = (4/ 2 )" 1 (p„ + 2p p ) with nuclcon number 
densities pi(i = n,p). It may be worth noting that the self-energy term a contains some divergences (see 
IV), and is reduced to a tre e — f~ 2 {^KpPp + ^KnPn) in the static limit of nucleons Jl4[ . 
Expanding the fluctuation field according to 

= \Zf E e 4(p - r+ "" r V4,5(n,p) (3.16) 
with the Matsubara frequency u> n = 2TmT, we find 

Z K f =[ II [^(n,p)][d&(n,p)]e i # / . (3.17) 

It is to b e noted that the zero mode (n = 0, p = 0) should be removed to avoid the double counting (see 



Eq. (2.16)). Here the effective action reads 



with the inverse thermal Green function, 

D eff = 



ai ( u l + {&- - Vi<)(u+ + P-K) + Pk sin 2 (0) 2(fl K + b)uj n \ , , 

P ^ -2(Ajc + &K m£ + (w_-/ijr)(« + +/i J r) ^ ' 

where 

J>± =±6+(p 2 + G 2 ) 1/2 (3.20) 

with G 2 = m^+b 2 and /1a" = Ma cos(0) with the effective mass, m^ 2 = cos(6)[m 2 K - — a]. It is to be noted that 
the effects of the condensate come in through the modifications of the chemical potential, fix — > cos (8), 
and the effective mass,mj^ = [m 2 K — a] — > cos{0)m*^ , except the additional term, y? K sin 2 (0). Then, the 
partition function can be simply represented as 

hxZ chiral =S C - Tr[G^{8G}] - l/2Tr'ln.D e// + Tr In Gh , (3.21) 

where the symbol Tr' means that the zero-mode should be removed. 

B. Dispersion relations for kaonic modes 

The excitation energy of the kaonic modes are given by the solutions w = iu n , which satisfy \D e ^\ = 
det(D eff ) = (see Appendix B): 

uj 4 - 2(ci + c 3 + 2c\)u? + 2cic 3 + cj = 0, (3.22) 
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where c\ = (w_ - + Mk), C2 

correspond to the if modes, 



+ 6 and C3 = \j1p? K sin 2 (6*). Then we have two solutions which 



Ei 



( Cl +c 3 + 24)±J(c 3 



2ci) 



4C1C 2 .. 



(3.23) 



Then 



b and C3 — > 0, with the dispersion relation in 
(lo± ± hk) 2 , which recover the previous 



In the limit (9) = 0, c\ — > (w_ 

the normal phase, w±(p) = ±6 + (p 2 + to J? + b 2 ) 1 ! 2 . 
results as they should do. 

In the condensed phase ((8) ^ 0) the mode corresponding to E- is the Goldstone mode as a consequence of 
the breakdown of V-spin symmerty; we can choose such form for the condensate as (K ± ) = f(8) / v / 2-exp(±z7) 
with an arbitrary 7. However, once a phase is chosen, the kaon-condensed state is a symmetry-broken phase 
with respect to the rotation around the third axis in V^-spin space by the U(l) operator, Uy(y) — cxpiiV^iy) 
with arbitrary while the effective Lagrangian is still inavariant under the transformation, U — > UyUUy. 
Thus we observe a spontaneous symmetry breaking (SSB) there. 

It is easy to show E_ ~ for p = 0. In the classical (tree-level) approximation Ci(p = 0) = (see 
Appendix B), which directly means E_(p = 0) = 0; actually we can expand E_, 



E 2 



_C3_ 2 

2C 2P 



P 

AC 2 



(3.24) 



for small momentum, |p| <C C. This dispersion relation may remind us of the Bogoliubov spectrum for the 
phonon mode in superfluidity or Bose-Einstein condensation |3(| . When we take into account the fluctuation 
effects (kaon loops) in the thermodynamic potential, their order of magnitude is O( h). A ccordingly, this 
relation should be modified in 0(K), because any kaon loop is not included in Eq. (3.23). ^ Hence we 
should find ci(p = 0) = + 0(h), which means that once fluctuations are included in the thermodynamic 
potential, the soft mode loses the Goldstone-boson nature. This situation is inevitable in the perturbative 
treatment [p7f : our thermodynamic potential takes into account the one-loop diagrams, whereas the self- 
energy diagram for kaons never include any kaon loop. However, we may expect its deviation to be small. 
Actually, we have found that the zero-point contribution from kaon loops is very small compared with 
tree-level contributions |T^,[l4|]. When we consider the thermal kaon loop, directly enters into the Bose- 
Einstein distribution function, /s(i?) = [exp((3E) — l] -1 , and it should diverge at p = 0. The other is the 
massive mode, E + 3> lOOMeV, and we may discard it for temperature we are interested in (T < lOOMeV). 
Thus only the thermal K~ loops play an important role due to this property. 

Moreover, if we can safely neglect the C3 term, we get a simple expression for E±, 



^4 



u)±±fi K = vV + C 2 ± (6 + Ak). 



(3.25) 



There are sever al reasons to support the pertinence of these formulae, especially for E-. First, we can see 
from Eq. ( |3.24| ) that in th e high momentum limit, |pl — ^sin(#), the difference of £L between 

Eq. (3.23) and Eq. ( p. 25 ) converges to be trifling. Q On the other hand, in the low momentum limit, 
|p| <C C, E- can be expanded as 



E^^ 
2C 



(3.26) 



within the tree approximation, which just corresponds to the second term in Eq. ( |3.24 ). Therefore, the 
approximated formula ( 3.25 ) may be well justified for the momentum, |p| 3> /i.R-sin(#), while it spoils the 
Goldstone-boson behavior (E- cx |p|) nea r p = 0. So we must carefully treat the small-p region, when we 
use the appro ximation. Comparing ( 3.23|) and ( 3.25 ) we can expect that the "nonrelativistic" dispersion 



relation ( 3.25| ) holds under the condition, v 2 = c 3 /2C 2 <C 1, namely the small velocity limit. Typically 
UK ~ O(100MeV) and C - 0(m K ) §, so that we find v 2 - 0(l(r 2 ) < 1. 

At finite temperature we can find a more meaningful criterion for the pertinence of Eq. (3.25). Consider 
a typical integral at finite temperature, 



11 On the other hand, the nucleon loops can be treated consistently (see pv| ). 

12 This property is very important for the renormalization of kaon loops, where ultraviolet divergences become 
relevant, see IV 
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d 3 p 



f B (E-(p)), 



which is the particle number of K fluctuations. Using Eq.( 3.25] ) we find 

(CT) 3 / 2 



(2tt) 3 / 2 



C(3/2) 



(3.27) 



(3.28) 



for T < C. However, as mentioned above Eq. ( 3.25 ) is no longer good for |p| < [Ik sin(#). We can estimate 
the integral for this region 



7T 2 /3 



(3.29) 



Therefore, the approximation ( 3.25| ) is verified, if Is /I <C 1 , which requires the following condition for C3 or 
Hk sin(0) as the function of density and temperature, 



C3 « S^C 2 (3/2) or m sin(0) « i(27rCT) 1 / 2 C (3/2). 
16 4 



(3.30) 



We can expect Eq. (3.30) for relevant densities and temperatures ]14fl , except very low temperature, where 
thermal effects should be trivially unimportant. 

Secondly, we can also expect the smallness of C3 qualitatively: in the limit, (8) — > or \lk — » 0, C3 gives 
no contribution. We also know that (9) and hk are inversely proportional to each other 0. Hence we might 
expect that the C3 term becomes relatively small. 



C. Thermodynamic relations 



The effective thermodynamic potential £l c hirai = —Tin Z c hi ra i reads from Eq. (3.21), 

^chiral — + ^sc + + ^Ni 

where fi c is the classical kaon contribution, 

He = V[-fm K (cos(e) - 1) - 1/2 • ii\f sin 2 (6)}, 



(3.31) 



(3.32) 



£l sc the subtraction term to avoid the double-counting of the interaction terms within the Hartree approxi- 
mation, 



Q sc = TTrlGff 1 {6G)] = V [acos(9)n K + bn 2 K ] . 
The kaon contribution, = — T\nZ e ^ , can be written as, 

where the zero-point energy (ZP) contribution fl K p and the thermal one f2^ are given as follows; 

n Z K P = V f -^kl[E + { V )+E_{ V )] 



and 



= TV 



d 3 p 



(2tt) 3 2' 



ln(l-e-^+(p))(l-e-^-(p)). 



The thermodynamic potential for nuclcons, iljv = —TlnZj^, can be also written as follows: 

Q N = -TTtIuGh = fii 



) zp + nth 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



(3.37) 



where ilf/* denotes the ZP contribution, 
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Q ZP 



-2V 



d 3 p 
(2^)3 



(3.38) 



and the thermal contribution, 



-2TVY^ 



n.p 



d 3 p 



ln(l + e -/3(^-w)) + i n (i + e -/»ta+«)) 



(3.39) 



Here, ej(e^), i ~ p,n are the single particle energies of the nucleons (anti-nucleons) , 

£ p,p = T Jp n K T - cos(0)) + £*, 
e„,» = Tj^k T ^(1 - cos<0» + E*, 



(3.40) 



with the kinetic energy, E* = \/p 2 + M* 2 . 
Using the thermodynamic relations, 



on 



chiral 



dT 



Qi 



on 



chiral 



E, 



chiral 



chiral 



■TS, 



chiral 



(3.41) 



we can find charge, entropy and internal energy. We easily se e tha t these quantities become too complicated 
to be tractable if we use the exact dispersion relation in Eq. (3.23) ; e.g. consider the kaonic charge given by 
way of the relation, Qk — ~dfl c hirai/dfiK, which results in a complicated form. Besides this, the quantities 
,n l K or n 2 Kl which can be written by the use of the thermal Green's functions, has a complicated form as 
well (see Appendix B). Hence a useful approximation is desir able. We have already discussed the pertinence 
of the approximation for the dispersion relation in Eq. (3.25) and we can see that this approximation allows 
us to write down the thermodynamic quantities in clear forms. For the kaonic charge it immediately gives 



)k/V = m f sin 2 (#) + cos(0)n^ + (1 + x) sin 2 ((0)/2) /9i 



where nx is the number density of thermal kaons, 

d 3 p 



n K 



(2tt) s 



[f B (E_(p))-f B (E + (p))} 



(3.42) 



(3.43) 



The resultant form ( |3.42 ) suggests that the charge of kaonic modes is screened by the condensate; in fact, 
we can see t hat t heir effective charge are given by ±e e // with e e // = cos(#)e by the use of the effective 
Lagrangian ( |3 . 1 5[ ) . Accordingly the modified chemical potential fin — cos{9)^lk gets the meaning of the 
chemical potential for the number of thermally excited kaons. 



IV. EFFECTS OF ZERO-POINT FLUCTUATIONS 



A. Renormalization 



In our formalism the thermodynamic potential ( |3.31 ), as well as the quantities n^, a, is constructed by 



many one-loop bubbles, which include intrinsic divergences. Hence we must renormalize them to extract 
finite contributions. Since the building units of these quantities are graphically one-loop bubbles of kaons 
and nucleons, it suffices to renormalize these diagrams pp| . Then all the diagrams with many loops have no 
more divergences. Since the temperature dependent part never induces new divergence, we renormalize the 
temperature independent part. Using the similar method proposed in ref. fl3|| , we can properly renormalize 
the thermodynamic potential. However, it is to be noted that we must carry out the renormalization program 
by including not only the kaonic contributions but also the nucleonic one in this case. The counterterms 
are introduced without spoiling the original symmetry structure of the Lagrangian; the symmetric terms 
are also invariant under £ and the symmetry-breaking terms should transform like Eq. (2.23). Since we 
have already known the stucture of the effective Lagrangian, given in Eqs. fl3.4|)-([3lj|), it is easy to find the 
proper counterterms. It is also to be noted that the vector- type contributions, SC induced by introducing 
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the chem ical potentials as well as the original Tomozawa- Weinberg term, never suffer from divergences (c.f. 
Eq. flC13| ), see also ref. |^9| ). Thereby it is almost clear the introduction o f che mical potentials never generates 
new divergences. In the following we use the approximated formulae ( 3.25| ) for the dispersion relations of 
kaonic modes for simplicity. 



First, we renormalize the self-energy terms of kaons and nucleons given in Eqs. (3.12) and (3.15), respec- 
tively. The kaon effective mass is given by 



1 K 



(4.1) 



The self-energy term a is given by the nucleon loop, er = /~ 2 Tr [GJ^jq'Ekn] = &tree + &ioop with the 
inverse propagator Gh (3.11). We only consider the loop correction <7i oop in the following. To t he one-loop 
order we don't care about the n x K that is given by the kaon loop in the effective mass M* (3.12). Then the 
temperature independent part 07 oop |t=o reads 

d 3 k M* 

1 / M* 2 

d 2 -M* 2 d + -M* 2 \n(^ T 



Vloop\T=Q — — 

i 



(4.2) 



with the cut-off (A) regularization, where d^s are the quadratically and logarithmically divergent integrals, 



d 2 = 
do = 



d 3 k 1 



A 2 



fe < A (2tt)3 2k 8n 2 ' 
d 3 k 11 
(2tt) 3 4fc 3 ~ 87T 2 



ln[^ 



fe<A 



(4.3) 



with an arbitrary momentum-scale k. To renormalize the effective mass we redefine the kaon mass and the 
KK interaction vertex. Then the counterterms should be given as 



o ct = 4 2jAf*E m .r 2 [a 2 - Af* 2 a ] 



(4.4) 



The meaning of the counterterms on is graphically clear (Fig. 4 ): ai renormalizes the kaon mass and ao 
the KK vertex at momentum p = 0. It is to be noted that the wavefunction renormalization is not needed 
to the one-loop order. 





(a) 



a 
(b) 



FIG. 4. Divergent diagrams (a) and the counterterms (b) to the one-loop order for the kaon propagator (dashed 
line). Solid line denotes the nucleon propagator without n} K and the cross (x) denotes the condensate (6). The black 
blob stands for the contribution from the counterterm. 



Then they are determined by the following conditions. 

d 2 



(Ploo P \T=0 + Cc 



= 0, 



T— @ct) \ i 



(4.5) 
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From Eqs. (pL2),(4.4)we find 



OL2 



M 2 
16tt 2 



a-o = d + 



1 


,„( 


'2k\ 


1" 








8^ 


k m) 


~ 2_ 



and consequently the finite contribution is given by 



T loop\T=0 



M* 2 ln 



M 2 



M 2 - M* 2 



(4.6) 



(4.7) 



Thus the self-energy term a is renormalized to be oy en , oyen = &tree + °i^£,|T=o- Note that the dispersion 
relation is modified by the replacement of a by a reni compared with the tree approximation where a is given 
by g — a t ree . Howe ver, we can see that the Goldstone-boson nature of is never spoiled by such 
procedure (see Eq. (4.27)). 



Next, we renormalize the nucleon effective mass, M*(i = n,p), which is given by Eq. ( 3.12 ). It includes 
the kaon loop correction n x K ^] , the temperature independent part of which is given by 



1 k \T=a 



d 3 p 



1 



(2tt) 3 E+ + E. 
d 2 ~ C 2 d 



(4.8) 



(see Eq. ( Cll )). Accordingly we introduce the counterterms, which renormalize the nucleon mass and the 
KN interaction vertex at momentum p = (see Fig. 5 ), 



Here C 2 is given by 



M l ct = /- 2 E«cos(#)[/3 2 - C 2 /3 ]. 



C 2 = (m 2 K - a) cos(0) + b 2 



(4.9) 
(4.10) 



but the nucleon loo p in a is irrelevant up to the one-loop order as w ell. It is to b e noted that the overall 
factor cos(0) in Eq. (|4.9|) stems from the transformation £ (see, Eqs. (3.4) and (|3.6|)). 



(a) 



Po 
lb) 



FIG. 5. Divergent diagrams (a) and the counterterms (b) to the one- loop order for the nucleon propagator. 
Notations are the same as in Fig. 4 . 



13 Note that another kaon-loop contribution n 2 K never includes divergence (see Eq. ( |C13[ )) because it originates from 
the vector-type KN interaction, while from the scalar-type interaction. 
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Then renormalization conditions are given as follows, 



{M*\ T=0 + M* ct )\ vac = M, 



9(0) 



: (M*\ T=0 + M* ct )\ va 



-A'i 



From Eqs. (3.12), 04.SD, (4.9) we find 



(h = d 2 - 



K 



16tt 2 ' 

and thereby the finite contribution to reads 

1 



1 




1" 




H-) 




8^ 


\m K ) 


~ 2 



1, ren. i 
n K \T=Q 



16tt 2 



m 2 K -C 2 - C 2 In 



c 2 



(4.11) 



(4.12) 



(4.13) 



It is to be noted that the structure of ( 4.12 ) is the same as (4.6). 

Once a and are thus renormalized, we can use the renormalized ones given by Eqs. (4.7) and (4.13) 
in the thermodynamic potential (3.31) instead of the unrenormalized ones. Although the thermodynamic 
potential (3.31) still includes the divergences through the ZP contributions in Eq. (3.35) and Eq. (3.38 ) , the 
total ZP contribution flzp = ^f" P + should be renormalized by the use of the counterterms aj and 0i 
without introduction of another counterterms (see Fig. 6 ). 



XX XX xxxx 




X X 

XX XX xxxx 



(b) 



(c) 



FIG. 6. Divergent diagrams to the one-loop order for the thermodynamic potential. Notations are the same as in 
Fig. 4. The diagrams in (a) are quadratically divergent and renormalized by the counterterms given in Figs. 4 and 5 
, while the diagrams given in (b) and (c), which are renormalized by redefining KK interaction vertex at momentum 
p = 0, give the quadratic and logarithmic divergences, respectively; the counterterms for (b) and (c) are already 
given in Figs. 4 and 5 . 



Thus we find the renormalized ZP contribution, £l^ e p, 



with the kaondoop contribution, 



64tt 2 



and the nucleondoop one, 

Wzp,n/ v = - ^2 rT-, 



32tt 2 



rjren r\ren , r^ren 



(C 2 - m 2 K )(m 2 K - 3C 2 ) + 2C 4 In 



(M* 2 - M 2 ){M 2 - 3M* 2 ) + 2M* 4 In 



M 



*2 



M 2 



(4.14) 



(4.15) 



(4.16) 



Here we have subtracted the vacuum energy which is quartically divergent, 



M 4 

ttvac/V = -3^4 - T^— o + 



I6ir 2 64tt 2 



(4.17) 



with the divergent integral, 



16 



/c<A 



d 3 k _ A 4 

(2tt)3 ~ 8^' 



(4.18) 



Then we can easily verify the relations, 



d(n r zp, K /v) 

dC 2 



K 



\T=0, 



2 S 



<9M* 



(4.19) 



It is to be noted that the magnitude of ZP contributions ( 4.15 ) and ( 4.16 ) are bounded from above, 

\W£p,kI V \ < 12MeV • fm- 3 and \Wj$ iN /V\ < 620MeV • fm" 3 . (4.20) 

Comparing the both limit values, we can see that the effect of the zero-point fluctuation by nucleons should 
be more pronounced than that by kaons. 

In the infinite-mass limit for nucleons, M -> oo, there is only left the kaon contribution and Q r £p is 
reduced to a simple form [l4}j , Q z e p — > tl r £p K . ^ ma y ^ e w °rth noting the relations, fi% — c?2 = fi and 
(/9b — do) + (/?2 — d2)/m 2 K = fo in terms of /, given in ref. jl3| , which is grap hically understood in Fig. 6 



Finally, we give a remark about the use of the approximated formulae (3.25) for th e dispersion relations of 
t he ka onic modes. As has been already mentioned in III, the difference between ( 3.25| ) and the exact formulae 
( |3.23 ) converges to zero for high- momentum limit, so that the structure of the ultraviolet divergences is th e 
same to each other. Therefore, even if we use the exact formulae, we find the same counterterms in Eq. ( 4.12| ). 



B. Nucleon-loop contribution 



We investigate the effects of zero-point fluctuations in neutron-star matter at T = 0, since the negative-sea 
effect is dominant there. At T = the thermodynamic potential is reduced into a simple form, 

Qchiral\T=0 = flc + Qsc\t=0 + ^Jv|t=0 + ^zp, (4-21) 

where f2 sc |T=o = V[a ren n 1 ^ ren \T=o cos(0)]. The nucleon contribution flN\r=o can be written as 

^at|t=o — V [3/5e£(l -x) + 3/5e p F x - /i„(l - x) - p p x] p B 

- {2bp K + a tree )f(l - cos(0)) (4.22) 

with the baryon number density ps, the proton mixing ratio x = p p /pB, and Fermi energies for nucleons, 

e F — (37r 2 ps(l — x)) 2 ^ 3 /2M , e p F = (3n 2 pbx) 2 / 3 /2M , in the nonrelativistic approximation. Since we have 
already known that the kaon contribution £l r z e p K gives only a tiny effect |T^Jl4[| , we, hereafter, concentrate 
on the nucleon one Wzp n by simply dropping kaon loops. After the calculation we will confirm it as a 
consistency check. 

It is well-known that the nuclear symmetry energy plays an important role for the ground-state prop- 
erties of the condensed phase. Hence we take it into account to get a realistic result ]14[] . Since we have 
already included the kinetic energy for nucleons, it suffices to consider only the potential energy contribution. 
Following Prakash et al. |3l"|] we effectively introduce a symmetry energy contribution, 

Q sy mm = VpB ^ _ 2x f S pot {u), (4.23) 

with the relative density, u = p B /po- The coefficient S pot (u) is given as S pot (u) = (S -(2 2 / 3 -l)(3/5)4)F(u) 
with the constraint F(l) = 1, to reproduce the empirical symmetry energy So ~ 30MeV, where e° F is the 
Fermi energy at p Q , e° F = (3n 2 p /2) 2 ^ 3 /2M , and F(u) is a simulated function, for which we take here 
F(u) = u for simplicity. 

Then the total thermodynamic potential, we consider here, is given by further adding the leptonic one 
^|t=o> 

&tot = Qchiral\T=0 + + &i\t=0, (4.24) 

where £Ii\t=o can be written in the standard noninteracting form, 
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fii|r=o = V 



Ac 

12tt 2 



o{\m\ 



YYY 

^[(2t 2 + l)Wt 2 + 1 - ln(f + Vt 2 + 1)] 



8?r 



3tt 2 



, (4.25) 



with i= J pr K — m^/m^. Here the first term is due to electrons and the second term due to muons, and we 

have used the relation yn = lik f° r leptons. 

The parameters (9) , x and fix are determined by the equilibrium conditions and charge neutrality of the 
ground state ||J||l3|. First, we get 



f sm(0)(m 2 



K G tre 



2 m b - & co S (0)) + £1 = 

d(0) 



(4.26) 



from the condition dVL to t/d{9) = 0. Recalling the relation in Eq. ( 1.19| ), we find it is further written as 

f 2 sin(0)(mf - 2(j, K b - p? K cos(0» = 0, (4.27) 

with the effective kaon mass mj? = m 2 K — cr ren . Then it is to be noted that the critical density is unchanged 
even if the zero-point fluctuation S!^™j, is included, since a ren is reduced to <7t ree as (8) — > 0. It is also to 
be noted that Eq. ( 4.27 ) ensures the Goldstone-boson nature of K~ mode (see Eqs. ( |3.23| ), ( 3.25 )), even if 
the ZP contribution is included. Secondly, we get 



4^ ot (u)(l - 2x) + 1 C 2 ° s(g) {2{V Kn - T,jcp) - li k } + m = 0. 



(4.28) 



from the condition dfltot/dx — 0, which is equivalent with the chemical equilibrium condition among kaons 
and nucleons, Lt n — fi p — lik- Finally, charge neutrality demands the relation, 



Qk + Qi = Q P (= xVp B ), 

which implies d£l to tl d^K — by definition, where the lepton charge is given by 

\ UK 



Qi/v 



/4 



3+3 



and the kaonic charge, 

Qk/V = n K f 2 sin 2 (fl) + {l + x) sm 2 ((9)/2)p B . 
Then the energy is given by way of the thermodynamic relation, 

Etot 



m,,t 



Qtot + f^iQi 
n,p,K,l 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



Note that these equations are the same as the previous ones in the heavy-nucleon limit MM if we simply put 



L ZP,N 



In Table I and II we present the optimum values of lik, x, and (9) for given densities, in comparison with 
the classical values. We use the values, aim s = —67 MeV, a 2 m s = 134 MeV and a^m s = — 134(— 222) MeV 
in Table I (II). [] 



Recent lattice simulations suggest larger values for KN sigma term [B2l 
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TABLE I. 


Optimum 


values of the parameters 


for the tree 


case and the 


one-loop case. We 


use the value, 


a :i m s = -134MeV. 
















tree 






tree + loop 




u 


(0>[deg.] 


fiK [MeVJ 


X 


Wldeg.J 


f-LK [Mev] 


X 


4.19 


0.00 


200.0 


n i no 
U.193 


0.00 


200.8 


n i no 
0.193 


4.40 


18.85 


244.0 


0.231 


18.74 


244.2 


0.231 


4.80 


31.97 


221.1 


0.293 


31.49 


222.6 


0.290 


5.20 


40.85 


19( . ( 


U.343 


on m 
39.91 


2U1.D 


0.OO8 


5.60 


47.62 


174.4 


0.384 


46.26 


181.3 


0.376 


6.00 


52.92 


151.9 


0.417 


51.26 


161.9 


0.407 


6.40 


57.09 


130.6 


0.442 


55.25 


143.5 


0.431 


6.80 


60.36 


110.8 


0.462 


58.46 


126.0 


0.450 


7.20 


62.91 


92.6 


0.477 


61.05 


109.7 


0.465 



TABLE II. Optimum values of the parameters for a^iris = — 222MeV. 



u 


(0>[deg.] 


tree 
Hk [MeV] 


X 


W[deg.] 


tree + loop 
Hk [MeV] 


X 


3.08 


0.00 


218.8 


0.156 


0.00 


218.8 


0.156 


3.20 


18.28 


206.4 


0.198 


17.42 


207.9 


0.195 


3.60 


39.75 


160.0 


0.326 


33.60 


181.3 


0.286 


4.00 


53.71 


109.2 


0.425 


42.20 


161.8 


0.344 


4.40 


63.74 


59.1 


0.494 


48.33 


144.7 


0.386 


4.80 


71.05 


14.4 


0.538 


53.12 


128.5 


0.419 


5.20 


76.13 


-23.3 


0.565 


57.04 


112.3 


0.444 


5.60 


79.67 


-54.9 


0.581 


60.30 


95.8 


0.465 


6.00 


82.19 


-81.5 


0.591 


63.08 


79.2 


0.482 
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Then we can observe a sizable effect of the zero-point fluctuation by nucleons; it shifts the optimum 
values of x and (9) downward, while [Ik upward, comapared with the classical values. The effect becomes 
prominent as the value of the KN sigma term is increased, since the reduction rate of the nucleon effective 
mass becomes high. As mentioned above it never affects the critical density, while it suppresses the growth 
of the condensate. 

In Fig. 7 the energy difference between the condensed phase and normal one is dipicted as a function of 
density. We can see that the energy gain is also reduced along with the suppression of the condensate; the 
zero-point fluctuation of nucleons reduces the energy gain by shifting the optimum values of parameters, 
while the resultant ZP contribution is small in magnitude. We also estimate the magnitude of the kaon loops 
by the use of the values given in Table I, II to be tiny for a consistency check. 




FIG. 7. Energy differences per baryon, AE/Nb = {E to t — E to t({0) = 0))/Nb, and ZP contributions for 
a3m s = — 134MeV (thin lines) and 03171, = — 222MeV (thick lines). Solid lines show the results with the zero-point 
fluctuation by nucleon, while dotted lines those in the tree level. Dashed lines show the ZP contribution. Kaon ZP 
contributions are estimated by the use of optimum values given in Table I, II (dash-dotted lines). 



We notice that these consequences originate from the replacement of the KN interaction term atree in the 



tree approximation by a ren in the effective mass of kaons 



the zero-point fluctuation works to reduce 



the attractive effect. It may be worth noting that m*£ stays to be in positive value and [Ik also becomes 
positive, once the zero-point fluctuation is included. It may remind us of the relativistic effect by nucleons 
p3[ , which also gives rise to the same effect. Hence, if we further take into account the relativistic effect 
in our formalism, it may be milder than before. Dispesion relations are also checked to be fairly good by 
comparing the results given by Eqs. |T23| ) and ( pT25| ). 



V. SUMMARY AND CONCLUDING REMARKS 



In this paper we have presented a formalism to treat fluctuations around the condensate within the 
framework of chiral symmetry. Our approach is based on the group theoretical argument; Goldstone fields 
are regarded as coordinates on the chiral SU(3)l x SU(3)r/ SU(3)v manifold and each state corresponding 
to the condensed phase or the normal phase with fluctuations can be represented as a chiral rotated state 
from the vacuum, which corresponds to the fiducial point on the manifold. Hence the fluctuations in the 
vicinity of the condensate are introduced by way of the successive chiral transformations by £ and rj. This 
procedure can be regarded as the introduction of the local coordinates around the condensed point on the 
chiral manifold. Since the chiral transformation is the isometry, the geometry around the condensed point 
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is the same as in the vicinity of the vacuum, where we might use the flat curvature. Hence this method has 
the advantage of avoiding the Lee- Yang term. 

More interestingly, this method corresponds to the separation of the zero mode; since there should appear 
the Goldstone mode as a result of the symmetry breaking of V-spin beyond the critical density, we must 
treat the zero modes nonperturbatively even in the normal phase. In the infinite- volume limit the condensate 
only contributes to the thermodynamical potential among zero modes in the condensed phase. Therefore we 
have put U into the form given in Eq. (2.14) from the beginning. It is worth mentioning that this method 
might have a wider applicability for studying phase transitions, e.g. in finite volume, if system has some 
symmetry. 

We have considered the thermodynamic potential up to the one-loop diagrams of kaons and nucleons and 
the KN interactions have been treated self-consistently within the Hartree approximation; consequently 
it consists of infinite series of many one-loops. Then the dispersion relation for kaon excitation, which is 
a fundamental object to examine the thermal properties of the condensed phase, can be obtained. There 
appear two modes: Goldstone-like soft mode and very massive mode, which correspond to K~- and K + - 
mesonic excitations, respectively, in the condensed phase. Hence, the thermal loops due to the soft mode 
play an important role. We have seen that the form of the dispersion relation can be reduced to a simple 
one, if we can take an approximation, C3 = 0. It should be worth noting that once this approximation works, 
the expressions of other quantities like charge, entropy, internal energy become very simple. 

The thermodynamic potential as well as the self-energy terms for kaons and nucleons contains some diver- 
gences. We have shown that these divergences are properly renormalized by way of the cut-off regularization; 
counterterms are introduced to redefine the masses of kaons and nucleons and the KN and KK interaction 
vertices. Once the self-energy terms are renormalized, wc need not require more counterterms to renormalize 
the thermodynamic potential. 

We have discussed the effect of the zero-point fluctuation by nucleons. It provides a sizable effect for 
the optimum values of the parameters, so that the energy gain is reduced from the classical values. This is 
because the KN sigma term remarkably reduces the nucleon effective mass in the condensed phase, while 
the critical density is not affected. On the other hand the contribution by the kaon loops is still tiny. 

It may be interesting to study thermal effects in the condensed phase, which is closely related with the 
problem of protoneutron stars, especially the effects of thermal loops of n 1 I ^ ren ,nj < Q). In a subsequent 
paper flS]] we discuss the thermal properties of the kaon condensed phase on the basis of the formalism given 
in this paper. We shall see the phase diagram in the density-temperature plane or equation of state for the 
kaon-condensed phase. 
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APPENDIX A: SEPARATION OF ZERO MODES 



Here we reconstruct our formalism from the viewpoint of separation of the zero mode, following ref. [ pTsj ] . 
To make our argument clear we first consider the finite volume V = L 3 and put it infinity at the end. Since 
we are interested in the kaon condensation, we res trict ourselves to the kaon and nucleon degrees of freedom 
for simplicity. Considering the partition function (^) and integrating out the baryonic degrees of freedom, 
we get the partition function for U, 



J chiral 



= N' J [dU]exp[S eff {U)}, 



(Al) 



with the effective action S e ^ \U) — S(U) + Tr In G(U), where S(U) stands for the pure mesonic action and 
the second term stems from the integration over the baryon field. The capital letter for trace is used for the 
Dirac operator G(U) of infinite dimension over space-imaginary-time. Substituting U = exp(2i<j> / /) with 
4> = T a 4>a in Eq. (Al), the action has the following structure, 



with 



S K = 



dr 



d 3 x 



The matrix G reads 
with the free part 

and the interaction part 
6G ^ KN 



S eff {U) = S K {K) + TrlnG, 



-d^K+d^K- - (m 2 K - n 2 K )K+K- - fi K K+ d T K~ + 0(\K\ 4 ) 



G = G + 5G 



Go 



d_ 



1 3 + r 3 



K+ d T K~ - 2^ K K + R- 



j 2 " j 2 g 

For the purpose of the Hartree approximation, we define Gr as 

G H = G + {6G), 



0(\R\ 4 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



by taking the thermal average for the kaon field. Then the effective action of the nucleonic part can be 
written as 



TrlnG = Tr In Gh + Tr [G^dG] - Tr [G^G)] + 0(\R\ 4 ). 
Thus effective action S e ^ can be written as 

S £ M = S K + TrlnGtf + Tr [G^JG] - Tr[Gt~^{6G)] + 0(|^| 4 ) 
Substituting Eqs. ((A§) into Eq. @), we find 
S eff = TrlnG// - TriGj/iSG)] 



(A8) 



(A9) 



dr / d A x 



-d^K+d^R- - M^K+R- - (b + n K )R + d T K~ + Q{\R\ 



(A10) 



K 



IIkat, with the self-energy term, Hkn = 
— o — 2[i K b. Here a and b stand for a = /~ 2 Tr[7 £if atG^ 1 ] and b — /~ 2 Tr[(3 + t 3 )/8G^ 1 ], respectively. It 
is to be noted that the "mass" term should include the contributions from the meson-baryon interactions 
(Hkn) and the chemical potential besides the genuine mass of kaons, m 2 K . 
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To perfo rm the integral in Eq. (Al) we expand the Goldstone field in terms of periodic plane waves like 

Eq. £m> , 



y^^ ex p[ 27 ™( n ' r / L " 



•n 4 rT)]C. 



(All) 



Inserting Eq. ( All ) into Eq. (A10), we can see that the magnitude of 4>a is at most 0(T/ \J + L~ 2 ). 
Hence if T <C Mk, the Gaussian integral should give the dominant contribution to the partition function. 

However, if the phase transition occurs and the "mass" term M%- tends to be vanished at the critical 
point, the standard chiral perturbation fails, because there is no quadratic term for the zero-mode (n M = 0); 
in fact, the "mass" term takes negative values before the phase transition, while it is vanished at the critical 
density. Hence we need to reorder the perturbation series by treating the zero mode nonperturbatively to 
get a consistent formulation valid even at the critical point. 

Following Gasser and Leutwyler we treat the zero-modes as the collective variables [pi, 



U = vllfV 

with Uf = exp[2i(T4?74 + T^ij^)/ /] and the U(l) matrix v — exp(i#oT4). 
representing the zero mode and r\ a sums up the nonzero modes, 



(A12) 

v is the constant unitary matrix 

(A13) 



After integrating out with resp ect t o the nonzero modes, we find the integral over the zero mode (U(l) C 
SU{3)) . Substituting Eq. (|AlJ) into Eq. (|A|), we can see that the form of the effective action is the 
same with Eq. (E^?l) with 9q instead of (9). Hence the integral over the zero mode becomes 



; = / d9 exp[-ft C M ra i(9 )/T} 



(A14) 



where f2 c 

volume dependence from the thermodynamic potentials, we write Eq. (A14) as 



is given by Eq. (3.31) with replacing the thermal ave rage (9) by 9q. Separating the trivial 



d9 exp[-V/3Lo] 



(A15) 



with uj = Sl c hirai{9o) /V . In the large volume limit or the low temperature limit we can apply the steepest 
descent method to evaluate the integral in Eq. ( Al5| ): 



: exp[-Wfo(0 



2tt 
VPui" ' 



where u>" stands for d 2 oj/d9Q\e 0= e™ and the saddle point 6*™ is determined by the equation, 

du(0o)/dO o \eo=er = 0- 



(A16) 



(A17) 



This is nothing but an extremum condition for the thermodynamic potential with respect to the order 
parameter. Hence 0™ is equivalent with the thermal average of 9, 9™ = (9). Accordingly the zero-mode 
matrix v rn is reduced to C- In conclusion, when some degrees of freedom should condense, it is enough to 
separate them from the beginning in the form Eq. ( A12| ) in the limit, V — > oo. 

Here it is also interesting to observe that there is a logarithmic singularity in the thermodynamic potential 
at the critical point as long as V ^ oo. 



15 In the general case the zero modes cover the SU(3) manifold, and the partition function is given by the group 
integral over SU(3) §cj. 
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APPENDIX B: CLASSICAL APPROXIMATION - (MODEL-INDEPENDENT) 

In the classical approximation, <j) a = 0, the operator u becomes trivial, 

U f = u = 1. (Bl) 

Furthermore, replacing the bilinear operator of the nucleon field by its expectation value, e.g., i/jip — > {4>ip), 
then we get 

SC = ^tr{flt[( C t [ Tem , C] + C[r em , C f ]), B]} + ^tr{[T em , C 2 ][T em , C f2 ]} 

- MA-(cos(0) - 1)(1 + x)p B - 1/2 ■ fi 2 K f 2 sin 2 (0), (B2) 

which stems from the kinetic terms of Goldstone bosons, and thereby model-independent ||. 
The symmetry breaking part is also transformed, 

Csb{U k = C\B) = C SB (1,B) + (cos(0) - 1)(Z K (1,B)) 

= (cos<0) - l)(-fm 2 K + Z KpPp + ^KnPu) + C SB (l, B), (B3) 



according to the octet dominance hypothesis (see, Eq. (|2.28| )), where we have used the nonrelativistic ap- 
proximation for nucleons. Therefore the result is also model independent once we employ this hypothesis 
||. Thus we can recover the previous result: for the energy difference between kaon condensed phase and 
normal one, 

8E/V = E cl /V + ti K Q? l /V, (B4) 

where the "classical" energy E c i, which is nothing but the classical thermodynamic potential, is given by 

E cl /V = (cos<0) - + x)p B /2 - fm 2 K + H Kp P P + ^KnPn)] - 1/2 • p\f sin 2 (#), (B5) 

for T = 0. The charge of the condensed kaons, is given by 

Qci/V = sm 2 (^) + ma-(1 + x)p B (l - cos(0})/2. (B6) 

It is to be noted that this result has been obtained without recourse to the definite Lagrangian. 

The optimum values of the parameters, (9),x,/j,k, are determined by the extremum conditions for E c i, 

dE cl /d{6) = 0, dE cl /dx = 0, dE cl /dp K = 0. (B7) 

In particular, we find the field equation for the condensate from the first equation, 

fsm{e){m*K - 2\x K b - y? K cos(0» = (B8) 

with m*K = m 2 K — <7 tree , which implies ci(p = 0) = for (0) ^ 0. 
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APPENDIX C: KAON GREEN'S FUNCTIONS 



The Green's functions of kaons can be extracted from the matrix D e f'. Define the inverse- propagator 
matrix V- 1 = /3~ 2 D eff : = + c x + 2c 3 , = u>l + a, and T>^ = -©21 = - 2c 2^n- Then the 

propagators are given by 

(Mi)=Vn=Vm\V~ 1 \~ 1 ( C1 ) 

and 

(<f> 5 cb 5 )=V 22 =V n 1 \V- 1 \-\ (C2) 

with det(2? _1 ) = |X> 1 1 . Recalling the Dzyaloshinskii relation [f23| . 

2?(w„) = D K (iw„), (C3) 

with the retarded (real-time) Green's function Z? fl , we find the excitation energy of the kaon modes are given 
by Dff 1 = or det(V- x ) = 0. 
Similarly 

(04^5) = ©12 = -V^\V- x \-\ and (0 5 4 ) = ©21 = -V^V- 1 ]- 1 , (C4) 



1. Calculation of n^n^- ,nf<- 



Consider the following quantities; = (K + K ), = ((K + K — K + K )) and nf^- = + 2fm 1 K . 
Hence it is sufficient to evaluate mP K and n x K . By way of the kaon Green's functions they are represented as 



4 = (K+K-) = l/2r J E / (^3(^11+^2) 

d 3 p 



(2tt) 3 



(C5) 



and 



= ((K+K- - k + k~)) = /r 1 £ / tct 2 ^!^" 1 !" 1 ^ 1 . 



(C6) 



Generally it is rather complicated to evaluate these quantities, but we shall see that we can easily do them 
under the approximation C3 = 0. In this case, T>^ = X>^ , thereby the determinant I© -1 ! is factorized into 
{V^ + iD^ 2 l ){V^ - ffl^). Hence we find 



det(V~ 



— r + — r ) (^ii 1 ) -1 



1 



1 



©n + ^12 ©u - *©12 



Note that T> x l ± iT> 12 x — [oj n ± iE^)(uj n T iE + ), then n° K , n 1 ^ can be written as 



^4 



E++E- \ L0l + £ 2 ^2 + £ 2 



(C7) 



(C8) 



and 



<i 3 p 



E++E- \u* + E' 2 _ ul + E'i 



(C9) 
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Using the relation 



E 



n— — oo 



X Z + n A 



• cotli7ra;, 



(CIO) 



we finally get 



' K 



d 3 p 



(2tt) 3 E+ + E_ 



1 



1 



(Gil) 



and 



d 3 p 



1 



(2tt)3 (£ + + £_)/2 



-(£_-£+)/2- 



£4 



(C12) 



The first terms in Eqs. ( Cll ) and (|C12j) mean the contributions from the zero-point fluctuation to be 
renormalized. Finally n 2 K — n° K + 2/i^n^ can be written as 



'K 



d 3 p 



(2tt)3 (E++E-)/2 



e pE + - 1 efi E - - 1 



(C13) 



It is to be noted that there is no contribution by the zero-point fluctuation in Eq. ( |C13| ) due to the vector 
nature of the KN interaction. 
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